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We have investigated, using effective mass approach (EMA), magnetic properties of a one- 
dimensional electron gas in graphene armchair ribbons when the electrons of occupy only the lowest 
conduction subband. We find that magnetic properties of the one-dimensional electron gas may 
depend sensitively on the width of the ribbon. For ribbon widths Lx = 3Mao, a critical point 
separates ferromagnetic and paramagnetic states while for Lx = (3M -I- l)ao paramagnetic state is 
stable (M is an integer and ao is the length of the unit cell). These width-dependent properties are 
a consequence of eigenstates that have a subtle width-dependent mixture of K and K' states, and 
can be understood by examining the wavefunction overlap that appears in the expression for the 
many-body exchange self-energy. Ferromagnetic and paramagnetic states may be used for spintronic 
purposes. 
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I. INTRODUCTION 

Generation of spin-polarized current is of significant 
importance both scientifically and technologically. A 
spin-polarized current would emerge naturally from fer- 
romagnetic materials. In order to achieve a spintronic de- 
vice, it is important to find non-magnetic materials where 
a spin-polarized current can be flowed without becoming 
depolarized. It has become possible to induce and detect 
spin polarization in non- magnetic semiconductors [lj-l3|l. 
In this paper we explore the possibility that graphene[J| 
field effect transistors [sj based on armchair nanoribbons 
can be used to generate spin-polarized currents. The 
conduction electrons in carbon-based materials can move 
very long distances without scattering due to their small 
spin-orbit coupling and low hyperfine interaction. Re- 
cently magnetic effects using graphene edges states of an 
armchair ribbon have been explored^. In our work we 
investigate bulk magnetic properties of a one-dimensional 
electron gas of armchair ribbons, where electrons in- 
teract via long-range Coulomb interaction. It is possi- 
ble that a ferromagnetic state may be stable in such a 
system as the investigation of one-dimensional Hubbard 
model with nearest neighbor electron-electron interaction 
suggests 0- 

Graphene armchair ribbons [sLfTol] have several special 
properties that are well suited for spintronic applications. 
According local density approximation (LDA) when the 
width of an armchair ribbon is = 3(M -I- l)an or 
= 3Mao a gap exists in the energy spectrum Q (a 
rather small gap exists when L^ = (3M -|- 2)ao, and this 
case will not be considered here). The other property is 
that boundary condition on the armchair edges admix K 
and K' valleys, and eigenstates are mixture of K and K' 
states forming one-dimensional subbands^] (this is in 
contrast to parabolic and cylindrical potentials [Tll - [l3| , 
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FIG. 1: (a) Electron density is such that the exchange self- 
energy is smaller than the Fermi energy and the electron gas 
is partially spin-polarized. We assume that, among conduc- 
tion subbands, only the lowest energy conduction subband 
is occupied with electrons. Dashed line indicates spin de- 
generate subband energy in the absence of electron-electron 
interactions. Spin-up and -down subbands are shown, (b) 
For a smaller electron density the exchange self-energy can 
be larger than the Fermi energy and the electron gas is fully 
spin-polarized. 



where the mixing is rather small). When the system 
is doped electrons occupy these subbands and a one- 
dimensional electron gas forms, see Fig[TJ The other 
unique property is a rather small value of the dielec- 
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trie constant (e ^ 1), which makes effects of Coulomb 
interaction effects strong. Many-body self-energies are 
thus one order of magnitude larger in comparison to the 
corresponding values of an electron gas of ordinary semi- 
conductors with e ^ 10. 
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large values of the width > 45A, where the LDA and 
EM A results agree approximately (see Sec. II). The width 
dependence of T,ex is shown in FiglJl Using these results 
we find that, for ribbon widths — 3Mao, a critical 
point separates ferromagnetic and paramagnetic states 
while paramagnetic state is stable for — (3Af + l)ao. 
These results are illustrated in FiglH This dependence 
on the width can be understood by examining the wave- 
function overlap that appears in the expression for the 
exchange self energy. The effect is a consequence of eigen- 
states that have a subtle width-dependent mixture of K 
and K' states. The large difference in the values of dielec- 
tric constants of graphene and ordinary semiconductors 
cannot explain these width-dependent magnetic proper- 
ties. 



FIG. 2: Excliange self energies of doped graphene for ribbon 
widths Lx = 3(M -|- l)oo (circles) and = 3Mao (squares), 
where ao = 2.46A and M is an integer. The effective mass 
approach is used. Here the dielectric constant is e = 3 and 
the dimensionless Fermi wavevector is kpao = 0.07. 






FIG. 3: (a) Bulk ferromagnetic state of one-dimensional elec- 
tron gas in a graphene armchair ribbon. We assume that only 
the conduction subband of a graphene armchair ribbon that 
is closest in energy to i5 = is occupied. The average spin 
value per carbon site is shown as vertical arrows (Only spins 
of electrons belonging to the lowest conduction subband are 
shown), (b) Phase diagram as a function of the Fermi energy 
Ef for the length of ribbon Lx = 3Mao. Ef is measured 
from the bottom of the conduction subband. F (P) stands for 
ferromagnetic (paramagnetic) state, (c) Phase diagram for 
Lx = (3M-I- l)ao. 



In this work we assume that a gap separates conduc- 
tion and valence subbands, and that, among conduc- 
tion subbands, only the lowest energy conduction sub- 
band is occupied with electrons. We calculate the many- 
body exchange self-energy of a spin-polarized one- 
dimensional electron gas of such a system for relatively 



II. EFFECTIVE MASS APPROXIMATION AND 
MODEL HAMILTONIAN 
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FIG. 4: Size of gap Eg of undoped graphene is shown. LDA 
results[l] for widths Lx = 3(M -I- l)ao (filled circles), Lx = 
3Mao (filled squares), and Lx = 3{M + 2)ao (filled triangles). 
EMA results for the same widths are open squares, circles, 
and triangles. 

Effective mass approach can describe numerous phys- 
ical properties of graphene approximately. It can be de- 
rived from tight binding method [T^. Fig U displays the 
value of the gap of armchair ribbons as a function of the 
width. Both LDA [9] and EMA results are shown for un- 
doped armchair ribbons. We see that the LDA values of 
the magnitude of the gap are in rough agreement with 
those of the EMA results, and the agreement becomes 
better for larger values of L^- However, there are sev- 
eral deviations. The EMA value of the gap is zero for 
Lx = (3M + 2)ao while that of LDA is small but non- 
zero. Also as Lx changes between 3{M + l)ao and SAfag 
the LDA value of the gap displays small oscillations while 
the EMA results do not. LDA includes electron-electron 
interaction effects and goes beyond nearest neighbor hop- 
ping while EMA does not. However, when L^ is larger 
than 45A and L^ = 3{M+l)ao or L^ = 3A/ao the EMA 
results are approximately correct. We will employ EMA 
to compute the approximate exchange self energy under 
these conditions. 
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We choose K and K' valleys as K 
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The ribbon is along y-axis. The enve- 



lope wavefunctions [ipA,ipB,] and [V'a'^b] '^^^ t)e com- 



bined into ^ 
Hamiltonian 
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V 



{lpA,->pB,-1pA^- 



tpg), which satisfies the 




Wavevector kx^y is measured from K and K' in the up- 
per and lower Hamiltonians, respectively. The hard wall 
boundary conditions 15] on each of the total A and B 
components of the wavefunction at the armchair edges 
X = and x — L'^ = + ao admix K and K' valleys. 
Note that envelope wavefunctions along the x-axis are 
g±ife„a; g^jjjj have the opposite wavevectors fc„ and —kn 
for K and K'. The wavefunction 5' of conduction sub- 
bands are 



^n{x,y,ky) 



I 



2JL 



^—ikji X 



(2) 



where 9{x) = 1 for < x < L'^ and otherwise. The 
wavevector of the n'th subband is fc„ = ^ — and 

" 3ao 

dk„,ky = Arctan(fcj^//c„). The eigenenergy is e„(fcy) = 
^ao^ky + k'^, where ^ — and t = 2.7eV. Subband 

energies near E = are plotted in FiglS] for two possi- 
bilities Lx = 3Mao and Lx = (3M -t- l)ao. In this figure 
the value of Lx is chosen so that in the first case fc„ > 
while in the second case kn < 0. 
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FIG. 5: Conduction and valence subbands with energies clos- 
est near E = 0. Solid lines are for Lx = 3Mao = 24ao. 
Dashed lines are for Lx = (3M -I- l)ao — 25ao. The numbers 
beside the curves indicate values of subband index n. 



III. EXCHANGE SELF-ENERGY 

We consider spin-polarized electrons in the lowest en- 
ergy conduction subband with the electron density no = 



— , where kp is the Fermi wavevector. In our model, we 
ignore the inter-subband mixing by Coulomb interaction, 
which is a standard approximation [l^. Then Hartree- 
Fock self-consistent eigenfunctions are plane wavesfll\ , 
given by Eq.Q, and no corrections of the kinetic energy 
are present. Note the Hartree self-energy will cancel with 
the potential of uniform positive background. 
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FIG. 6: Magnitude of dimensionless exchange self-energy 
A{Yf, Xn) as a function of Yf for X„ — n/3 (solid line) and 
—tt/3 (dashed line). Inset: magnitude of the dimensionless 
self-energy A{YF,Xn) as a function of Yf computed using 
Eq.Q (solid line) and Eq.© (dashed line) for X„ = tv/3. 



We denote these wavefunctions by jfcj,). At fcy = the 
exchange self-energy [l3| is 



S„,e.(0) = / dky{0,ky 



fc,,0 , 



(3) 



which can be written as T,n.ex{0) = —-^A{Yp, Xn) with 
the dimensionless exchange self-energy defined as 



A{YF,Xn) 

= ± r dx r (i 



(1 - cosX) 



(4) 



Here Y = kyL'^ 
dimensionless. 



Yp = kpL'^ and Xn — knL'^ are all 
In the case Lx — 3Mao we have Xn = 



7r(n — 2M — |), and the conduction subband with the 
lowest energy has the value Xn = 7r/3 corresponding to 
the value n = 2M+\. For Lx = (3M-f l)oo we have Xn = 
7r(n — 2M — |), and the conduction subband with the 
lowest energy has the value A"„ = —it/?i corresponding to 
n = 2M + 1. In Fig^A(Yp,Xn) is plotted as a function 
of Yp for Xn = 7r/3 and — 7r/3. For ^ 1 we can 
approximate 



A{Yp,Xn) 




(5) 
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where G is Meijer G- function [Ts!]. In Fig l6] this approxi- 
mate analytical result is compared with the exact numer- 
ical result of Eq.(|4]). As expected, the agreement between 
the two are excellent for small values of Yp . 

Exchange self energies of doped graphene for ribbon 
widths Lx = 3(A/ + l)ao and Lx = 3Mao are shown in 
Figl21 Here the dimensionless Fermi wavevector kpao — 
0.07. We notice that the self energy displays oscillations 
as Lx changes between 3(M-|- l)ao and SAfap. A similar 
effect was also observed in the LDA result for the gap in 
the undoped case (compare Figs H] and S]). 

Spontaneous spin splitting will occur in the lowest con- 
duction subband when |S„,ea:(0)| > Ep, see FigU] This 
condition is equivalent to A(yp, X„) > + Yp — 

\Xn\. At Xn — 7r/3 the inequality is satisfied when 
Yp < 11.309 for e = 1 and Yp < 3.22 for e = 3. The 
critical point between ferromagnetic and paramagnetic 
states is shown schematically in FiglHJb). At Xn = — 7r/3 
no critical point exists since A{Yp,Xn) is too small, see 
Figini In this case paramagnetic state is stable, as illus- 
trated in Figl3l[c). 

Also we must require that the exchange self en- 
ergy correction be smaller than the gap |S„^e2;(0)| < 
Eg = 27ao|fc„|. This condition is equivalent to 
^A(Ff,X„) < 2|X„| = f . We find at X„ - n/3 
that this condition is satisfied when Yp < 0.364 for e = 1 
and Yp < 2.439 for e = 3. 

The dependence on the width can be understood by 
examining the wavefunction overlap that appears in the 
expression for the exchange self energy. It is given by 
|V';^(a;,y, 0)V'n(a;, y, oc 1 cos6'fc„,fc„ for ri = r2. 

When fc„ < the quantity cos 6h h = , < 0, 

which makes the dimensionless self energy [s.{Yp,Xn) 
smaller, see the integrand in Eq.Q. This effect is a con- 
sequence of the presence of the phase factors e"**''" '"!' 
only in certain components of the eigenstate wavefunc- 
tions (see Eq.Q) and that the value of A;„ for the lowest 
conduction subband is negative when Lx = (3M -|- l)ao 
(see arguments below Eq.^). In other words it is a con- 
sequence of a subtle width-dependent mixture of K and 
K' states in the eigenstate wavefunctions. 

IV. SUMMARY AND DISCUSSIONS 

As the comparison with the LDA result shows, when 
Lx large and for Lx — 3(Af -I- l)ao and Lx = 3A/ao the 



result of EMA is approximately correct. Under these 
conditions we have employed EMA to compute the ap- 
proximate exchange self energy. We have shown that 
when only the lowest conduction subband of a graphene 
armchair ribbon is occupied magnetic properties of the 
one-dimensional electron gas may depend sensitively on 
the width of the ribbon. We find that, for ribbon widths 
Lx = 3Mao, a critical point separates ferromagnetic and 
paramagnetic states while paramagnetic state is stable 
for Lx = (3M -I- l)ao- This dependence is in sharp con- 
trast to one-dimensional electron gas of ordinary semi- 
conductors, and can be understood by examining the 
wavefunction overlap that appears in the expression for 
the exchange self energy of an armchair ribbon. It reflects 
the fact that eigenstate wavefunctions of an armchair rib- 
bon contain a non-trivial width-dependent mixture of K 
and K' states. The large difference in the value of dielec- 
tric constant between graphene and ordinary semicon- 
ductors alone cannot explain our width-dependent mag- 
netic properties. The magnitude of the exchange spin 

2 

splitting is of order ^^p-, and, depending on the size of 
the Fermi wavevector kp, it can vary in the range of 
10 — lOOmeV. As shown in FiglT]our work suggests that 
the degree of spin polarization in a graphene armchair 
ribbon may be controlled by changing the gate voltage. 

An estimation, using the usual mean field value, gives 
that Curie temperature ksTc is of the order of the ex- 
change self energy O.leV. However, quantum fluctuations 
may reduce this value. A spin-polarized CDW and a Lut- 
tinger liquid are possible true groundstates. It may be 
worthwhile to investigate these issues using improved ap- 
proximation schemes than EMA employed in this paper, 
such as spin density functional method and exact diag- 
onalization techniques. Our results suggest that exper- 
imental investigations of magnetism in one dimensional 
electron gas of an armchair ribbon may produce numer- 
ous interesting results. 



Acknowledgments 

We thank A.H. MacDonald for several useful sugges- 
tions. This research was supported by Basic Science Re- 
search Program through the National Research Founda- 
tion of Korea(NRF) funded by the Ministry of Education, 
Science and Technology (20I2R1A1A2001554). 



[1] B. V. Wees, Nature Phys. 3, 147 (2007). 

[2] S. P. Dash, S. Sharma, R. S. Patel, M. P. de Jong, and 

R. Jansen, Nature 462, 491 (2009). 
[3] S.-R. Eric Yang and A. H. MacDonald, Phys. Rev. B 67, 

155202 (2003). 

[4] A. K. Geim and A. H. MacDonald, Phys. Today, 60, 35 



(2007). 

[5] F. Schwierz, Nat. Nanotechnol. 5, 487 (2010). 

[6] H.-H. Lin, T. Hikihara, H,-T. Jeng, B.-L Huang, C.-Y. 
Mou, and X. Hu, Phys. Rev. B 79, 035405 (2009). In this 
work the edges states that are responsible for magnetism 
are partially filled. In our case these edge states are either 



5 



empty or completely filled and are not relevant. 
[7] S. Daul and R.M. Noack, Z. Phys.B 103, 293 (1997); 
R. Strack, D. VoUhardt, J. Low Temp. Phys. 99, 385 
(1995). In one-dimensional Hubbard model the ferro- 
magnetic state is unstable when only the on-site (short- 
range) interaction is used. However, Strack and VoUhaxdt 
demonstrate that the ferromagnetic state can become 
stable when the nearest neighbor Coulomb interaction 
is included. 

[8] A. H. Castro Ncto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

[9] Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 

Lett. 97, 216803 (2006). 
[10] L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 
[11] P. S. Park, S.C. Kim, and S.-R. Eric Yang, Phys. Rev. 

Lett. 108, 169701 (2012). 



[12] P.S. Park, S.C. Kim, and S.-R. Eric Yang, J. Phys.: Con- 
dens. Matter 22, 375302 (2010); P.S. Park, S.C. Kirn, and 
S.-R. Eric Yang, Phys. Rev. B 84, 085405 (2011). 

[13] S.C. Kim, and S.-R. Eric Yang, J. Phys.: Condens. Mat- 
ter 24, 195301 (2012). 

[14] T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

[15] C. Tang, W. Yan, Y. Zheng, G. Li, and L. Li, Nanotech- 
nology 19, 435401 (2008); S.C. Kim, P.S. Park, and S.-R. 
Eric Yang, Phys. Rev. B 81, 085432 (2010). 

[16] H. Haug and S. W. Koch, Quantum theory of the op- 
tical and electronic properties of semiconductors (World 
Scientific, Singapore, 1990). 

[17] N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Brooks Cole, 1976). 

[18] H. Bateman and A. Erdlyi, Higher Transcendental Func- 
tions, Vol. L (McGrawHill, New York, 1953). 



